Quantitative measurements of the thermal resistance of Andreev interferometers 
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Using a local thermometry technique, we have been able to quantitatively measure the thermal 
resistance 7?^ of diffusive Andreev interferometers. We find that is strongly enhanced from 
its normal state value at low temperatures, and behaves non-linearly as a function of the thermal 
current through the sample. We also find that the oscillates as a function of magnetic flux 
with a fundamental period corresponding to one flux quantum "l>o = h/2e, demonstrating the phase 
coherent nature of thermal transport in these devices. The magnitude of is larger than predicted 
by recent numerical simulations. 

PACS numbers: 74.25.Fy, 74.45. +c, 73.23.-b 
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A diffusive normal metal (N) in proximity to a super- 
conductor (S) in a mesoscopic hybrid device also acquires 
superconducting properties by the process of Andreev re- 
flection 0: at temperatures well below the gap of the 
superconductor, fc^T << A, an electron in the normal 
metal cannot be transmitted through the NS interface, 
but is reflected as a coherent hole with the simultaneous 
generation of a Cooper pair in the superconductor. The 
electrical transport properties of such proximity-coupled 
systems have been extensively investigated both experi- 
mentally and theoretically in the last decade 
More recently, the thermal transport properties have at- 
tracted much theoretical interest 0, 0, IS Q ; following 
measurements of the phase -dep endent thermopower of 
Andreev interferometers [lllll^illjl- Inelastic-scattering 
dominated normal-metal systems, the ratio of the elec- 
trical to the thermal resistance is proportional to the 
temperature, the so-called Wiedemann-Franz (WF) law 
However, theoretical studies indicate that the ther- 
mal resistance of a normal metal in the proximity regime 
is strongly enhanced 0, 1^, 13- The WF law, which is 
widely valid in a disordered metal system, is no longer 
correct for proximity-coupled systems. This topic was 
flrst explored experimentally by Dikin et al. 15]. How- 
ever, the sample in that experiment had two supercon- 
ducting elements directly in the path of the thermal cur- 
rent, so that the reduced thermal conductance may have 
been due to the well-known suppression of thermal con- 
ductance in a conventional superconductor 0, lisf rather 
than a proximity effect phenomenon. 

In this Letter, we report measurements of the tem- 
perature and magnetic fleld dependence of the thermal 
resistance i?-^ of Andreev interferometers (consisting of 
a hybrid loop with one superconducting arm and one 
normal-metal arm) without superconductors in the ther- 
mal current path, and therefore in the true proximity 
regime. As predicted by theory 011, 0|j we find that 
W of a proximity-coupled normal metal is enhanced by 
as much as one order of magnitude from its normal state 
value at low temperatures. also oscillates as a func- 
tion of the external magnetic flux, with a fundamental 



period of one superconducting flux quantum $o — h/2e. 
These oscillations demonstrate the phase coherent nature 
of the thermal current in this system. Furthermore, 
is a strongly nonlinear function of the thermal current l'^ 
even at thermal currents as small as a few femtowatts. 

Experimentally, the thermal resistance is deflned as 
the ratio of the temperature differential AT across the 
sample to the thermal current /'^ through the sample, 
i.e., R^ = AT//-^, under the condition that no electrical 
current flows through the sample (/ = 0). Hence, in or- 
der to obtain RF , our devices should be designed so that 
and AT can be measured directly. Figure 1(a) shows 
a scanning electron microscopy (SEM) image of one of 
the devices we measured. There are two layers of metal 
on this device. The first layer is a 50 nm thick Au film, 
followed by a 100 nm thick Al film deposited on top of 
the Au in a second level of lithography, after an in situ 
oxygen plasma etch was used to clean the Au surfaces in 
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FIG. 1: (a) Scanning electron microscope (SEM) image of 
one of the devices. The device consists of three parts: (1) On 
the left, a metallic heater line with a thermometer. (2) On 
the right, a large normal- metal pad with another thermome- 
ter. (3) In the middle, a 'house' Andreev interferometer, (b) 
Schematic of Andreev interferometers we measured with two 
different geometries: (left) 'house' and (right) 'parallelogram.' 
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order to obtain good NS interfaces. The device consists of 
three parts: (1) In the center is the sample itself, which in 
this case is an Andreev interferometer in the 'house' con- 
figuration (in the terminology introduced in Ref. 
The Andreev interferometer includes a one-dimensional 
(ID) normal- metal (Au) wire and a superconducting (Al) 
loop. (2) On the left is a heater, which is essentially a 
0.68-/im-wide Au line connected electrically to the sam- 
ple. By passing a direct (dc) current through this heater, 
one can raise the electron temperature at one end of the 
sample. (3) On the right is a large, normal-metal pad 
that serves as the cold end of the sample. Attached to 
both the heater and the cold normal-metal pad are ther- 
mometers that measure the local electronic temperature. 
The normal-metal heater is approximately 25 /zm long, 
much longer than the inelastic electron scattering length, 
so that one can define a effective local electronic temper- 
ature in the middle of the heater that can be measured by 
the thermometer on the left, which we denote the 'hot' 
thermometer (the thermometer attached to the other end 
of the sample is the 'cold' thermometer). The design and 
operation of these thermometers have been described in 
detail elsewhere and will not be discussed here. Con- 
nections to both end of the heater are made through su- 
perconducting Al contacts to reduce the heat flow. Two 
additional probes at each end allow us to measure the to- 
tal four-terminal differential resistance of the heater for 
any value of the dc current through it, thereby deter- 
mining the power generated in the heater. In total, five 
devices were measured, in two different configurations, 
the 'house' and 'parallelogram' geometries shown in Fig. 
1 (b) . Data for three of these devices are discussed in this 
paper. 

To determine the thermal resistance, the electrical re- 
sistances of both the hot and the cold thermometer are 
measured as a function of the temperature of the dilution 
refrigerator mixing chamber, with no dc current through 
the heater. The resistances of both thermometers are 
then measured as a function of the dc current Ih through 
the heater, with the base temperature of the dilution 
refrigerator fixed. The two measurements are then cross- 
correlated to obtain the effective electron temperature as 
measured by the hot and cold thermometers as a function 
oi Ih] a simultaneous measurement of the differential re- 
sistance of the heater determines the equivalent power Ph 
generated in the heater. The difference between the tem- 
peratures measured by the two thermometers gives AT. 
Since all connections to the heater except the connection 
to the Andreev interferometer are made via supercon- 
ducting contacts whose thermal conductance is negligible 
at temperatures far below the transition temperature, the 
power generated in the heater can only flow out through 
the Andreev interferometer, or through the substrate. At 
temperatures below about 200 mK, the electron-phonon 
coupling in the normal metal is very weak, so that ther- 
mal leakage to the substrate can be ignored |0| . In addi- 



tion, while our earlier devices (such as the parallelogram 
interferometer discussed in this paper) were fabricated 
on oxidized Si substrates, more recent devices (such as 
the house interferometer shown in Fig. 1(a)) are fabri- 
cated on 50 nm thick Si3N4 substrates. In these samples, 
the thermal leakage to the substrate is even smaller, so 
that in principle, the valid range of measurement can be 
extended to higher temperatures. In practice, the tem- 
perature range is restricted by the sensitivity of the ther- 




FIG. 2: Thermal resistance ii as a function of the heater 
power Ph at six different base temperatures for (a) the 'house' 
and (b) the 'parallelogram' interferometers, (c) vs. Ph at 
small values of Ph for the 'house' geometry. The solid lines 
are fits to the functional form oc ^Jl/Ph- For the 'house' 
interferometer, the distances from either side of the sample 
and from the NS interfaces to the center node are (referring 
to Fig. 1(b)) L = 1.55 nm and L' — 0.29 fim respectively. 
For the 'parallelogram,' L — 1.19 jim, L' = 0.24 fim, and 
L" = 0.76 nm. 
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mometers. Since all the power generated in the heater 
flows through the Andreev interferometer, the thermal 
current through the sample is simply given by the 
measured value of Ph, and the thermal resistance is 
given by AT/Ph- This measurement is repeated at dif- 
ferent values of the base temperature of the refrigerator 
to obtain i?^ as a function of temperature. 

Figure 2(a) and (b) show the thermal resistance 
of a 'house' and a 'parallelogram' interferometers at six 
different temperatures as a function of the power of the 
heater Ph- As shown in Fig. 2, the qualitative behavior 
of R^ as a function of Ph is similar for the two geome- 
tries. In order to experimentally define a thermal resis- 
tance in the linear response regime, i?^ should approach 
a limiting value as — > 0. However, R^ continues to 
change as a function of Ph (or equivalently, l'^) even at 
a heater power of a few fenitowatts. At low values of 
Ph, R^ oc ^Jl/Ph, as shown in Fig. 2(c). This power 
law dependence of R?^ on Ph is valid for all samples and 
at all temperatures measured. Theoretically, of course, 
a linear response R^ can be defined; however, numerical 
simulations for these devices based on the quasiclassical 
theory also find that R^ oc yJlJPh for intermediate 
values of Ph- These simulations show that the linear re- 
sponse regime is approached for Ph ^10 fW, comparable 
to the minimum values used in the experiments. 

Given the nonlinear behavior of R^ on Ph, we shall 
experimentally define the linear response thermal resis- 
tance as the value of R^ at the lowest heater power mea- 
sured. This corresponds to Ph = 3.7 fW for the 'house' 
interferometer and Ph = 31 fW for the 'parallelogram' 
interferometer of Fig. 2 respectively. The solid symbols 




0.06 0.12 
T(K) 



0.18 0.03 0.06 0.09 0.12 0.15 
T(K) 



FIG. 3: The solid symbols are the thermal resistance of (a) 
the 'house' interferometer and (b) the 'parallelogram' inter- 
ferometer as a function of the mixing chamber temperature, 
determined from the data of Fig. 2 at the lowest heater power 
applied in the experiments. The dotted lines are guides to the 
eye. The solid lines represent the thermal resistance of equiv- 
alent normal-metal wires, estimated using the Wiedemann- 
Franz law, and the measured normal-state electrical resistance 
of the interferometers. The dashed lines represent theoretical 
calculations of the thermal resistance of the Andreev interfer- 
ometers, using the experimental parameters for the samples, 
as described in the text. 



in Figs. 3(a) and 3(b) show the resulting thermal resis- 
tance as a function of temperature for the 'house' and 
'parallelogram' interferometers. For comparison, we also 
show the expected thermal resistance for an equivalent 
normal-metal sample, calculated using the Wiedemann- 
Franz law from the measured normal-state resistance of 
the wire, using the textbook value for the Lorenz number. 
For both samples, R^ increases rapidly with decreasing 
temperature. For the 'house' thermometer, R^ is larger 
than the thermal resistance for an equivalent normal sys- 
tem by almost an order of magnitude at the lowest tem- 
perature. For the 'parallelogram' interferometer, at first 
sight, the increase is not as large; however, it must be 
noted that R^ for this sample was inferred at a heater 
power of 31 fW. If we extrapolate R^ to a value of 3.7 
fW as for the 'house' interferometer, we obtain values of 
R^ of the same order of magnitude as in Fig. 3(a). 

Depending on the dimensions of the sample, the elec- 
trical resistance of proximity coupled normal-metals can 
show 'reentrant' or non-monotonic behavior, where the 
resistance first decreases below the transition of the su- 
perconductor, but reaches a minimum and then starts in- 
creasing as the temperature is decreased further Q . The 
non-monotonic behavior is associated with the competi- 
tion between two effects, a decrease in the resistance due 
to pair correlations induced in the normal metal, and an 
increase in the resistance due to a decrease in the density 
of states N{E) near the Fermi energy Ep. At any finite 
temperature, the effect of pair correlations is greater, but 
at T = there is expected to be an exact cancellation 
for systems without interactions, resulting in the system 
regaining its normal-state resistance at T = 0. Theo- 
retically, the thermal resistance of the 'house' interfer- 
ometer is expected to be influenced only by the decrease 
in N{E), and hence increases monotonically as the tem- 
perature decreases. Unlike a superconductor, however, 
N{E) does not go to zero as T ^ 0, but saturates at 
a value that depends on the dimensions of the sample 
and the transparency of the NS interfaces. As T ^ 0, 
N{E) is small but finite; the quasiparticles occupying 
the levels in the pseudogap contribute to the thermal 
conductance, leading to an enhanced thermal resistance, 
but one that still varies inversely with T according to 
the Wiedemann-Franz law. Numerically, i?"^ approaches 
this limiting behavior below a temperature correspond- 
ing to approximately OAEc/ks, where Ec = TlD/Lq is 
the correlation energy 8] . Here D is the electronic diffu- 
sion coefficient in the normal metal, and Lq is the length 
from one end of the interferometer to one of the NS in- 
terfaces. In the 'house' interferometer of Figs 2 and 3, 
for which the diffusion constant D = 208 cm^/sec, and 
Lq = 1.84 /j,m, O.lEc/kB ^ 4.7 mK, below the tempera- 
ture range of the experiment. Equivalent parameters for 
the 'parallelogram' interferometer are D = 127 cm^/sec 
and Lq = 1.43 /im , giving a similar saturation tempera- 
ture of O.lEc/kB ^ 4.7 niK. The dashed lines in Fig. 3(a) 
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FIG. 4: The solid symbols are the thermal resistance of a 
'house' interferometer at different magnetic fields, measured 
at T = 40 mK. The solid line is a guide to the eye. The 
dashed line is the resistance of the same sample as a function 
of the magnetic field at T = 400 mK. 



and 3(b) show the result of the numerical simulations of 
the thermal resistance of the 'house' and 'parallelogram' 
interferometers respectively, with the parameters given 
above, and assuming perfectly transparent NS interfaces 
The theoretical predictions show significant devia- 
tions from the normal state thermal resistance only at 
temperatures below 20-30 mK, while the experimental 
i?"^ is already larger than the normal state thermal re- 
sistance at temperatures about an order of magnitude 
higher. We believe this deviation might be due to an 
intrinsic mesoscopic effect not restricted to NS devices 
associated with the long length scales required to equili- 
brate the energy of the quasiparticles in mesoscopic sys- 
tems. 

The phase dependent nature of transport in mesoscopic 
NS structures means that the thermal resistance of an 
Andreev interferometer should oscillate periodically with 
an externally applied magnetic field, with a fundametal 
period corresponding to a flux h/2e through the area of 
the interferometer loop 0, Detailed simulations in 
Ref. show that the 'house' interferometer has a larger 
oscillation amplitude of the thermal resistance compared 
with the 'parallelogram' interferometer, hence we focused 
on the former in measuring oscillations of . Experi- 
mentally, measuring the oscillations in turns out to 
be a difficult proposition. Ideally, we would like to bias 
the sample at a specific heater power and look at the vari- 
ations in the temperature measured by the local electron 
thermometers as we sweep the magnetic field. However, 
at a heater power of 3 fW, the maximum variation in the 
temperature difference would be on the order of 0.5 mK 
(as can be seen from Fig. 2(a)), below our experimental 
resolution. One could apply a larger amount of power 
to the heater; however, the maximum variation in R"^ 
drops dramatically at larger values of Ph- Consequently, 
in order to measure the variation of R^ with external 



magnetic field B, we use the same technique as for the 
temperature dependence: at a fixed external magnetic 
field B, we measure the thermal resistance as a function 
of heater power P/j, and take the value at the lowest mea- 
sured heater power as the value of RF at that value of 
magnetic field B. Figure 4 shows the resulting R^ for a 
third Andreev interferometer device in the 'house' con- 
figuration, measured at T=40 mK, at Ph= 2.1 fW. For 
comparison, we also show oscillations of the resistance 
of the interferometer measured at r=400 mK. It can be 
seen that R^ oscillates periodically as a function of B, 
with a period corresponding to one superconducting flux 
quantum $o = h/2e = 0.092 mT through the interferom- 
eter loop. Both the electrical resistance R and RF have 
the same symmetry with respect to B. Since R is known 
to be symmetric with respect to -B, this means that RF is 
also symmetric with respect to B, as predicted by theory 
(the offset seen in the data is due to the remanent field 
in the external superconducting solenoid). 

In summary, we have measured the thermal resistance 
of Andreev interferometers in two different geometries. 
We find that the measured thermal resistance of all the 
samples is enhanced at low temperatures and deviates 
from the values estimated from the WF law for equiv- 
alent normal-metal systems. In addition, the measured 
thermal resistance shows strong nonlinear behavior with 
respect to the thermal current . At small values of J-^, 
R^ (X y^l/I'^. Finally, we have observed that the ther- 
mal conductance oscillates periodically as a function of 
the applied magnetic flux with a fundamental period cor- 
responding to <I>o = h/2e. Our results are qualitatively 
consistent with recent numerical simulations 0. 

This work is supported by the NSF through grant num- 
ber DMR-0201530. 
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